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ABSTRACT 


This is the second in a series of reports on the use of tracking information from NASA's 
Syncom El satellite to define the longitude-dependent gravity field of the earth. This tracking 
information shows that Syncom n drifted freely in the gravity fields of the earth, sun and moon 
from 25 April to 3 July 1964. In early May, the ascending Equator crossing of the orbit, as 
determined from radar range and range-rate tracking at the Goddard Space Flight Center, was 
at 121° West, with a westward drift rate of 0.81 degrees/day. During the following two-month 
period, this drift rate decreased steadily until late in June when, with the ascending Equator 
crossing at 161° West, the westward drift rate was 0.75 degrees/day. This decelerating drift 
rate is entirely consistent with the magnitude and phase angle of the equatorial ellipticity of the 
earth, as previously determined by this author. These ellipticity constants were responsible for 
the accelerated drift of this satellite in the longitudes from 55° to 64° West. By utilizing the 
first (energy) integral of the drift motion of a 2 4- hour satellite in a second-order longitude- 
dependent earth gravity field, the drift of Syncom n over the Central Pacific during the two 
months in question was seen to be sensitive to the following parameters of the earth's equatorial 
ellipticity. 

J 22 = -(1.71 + 0.22) x 1CT 6 , 


which corresponds to a difference in the major and minor equatorial axes of the earth of 
a 0 - b 0 = 214 ± 28 feet, and 


A. ^ 2 - -( 17. 1 ± 4.9 ) degrees , 


which gives the longitude of the major equatorial axis with respect to Greenwich. These new 
results confirm the postulate (made in an earlier study) that higher order longitude- dependent 
earth gravity has small influence on the long term drift of the high- altitude 24- hour 
satellite. The consistency of the new drift data indicates that the earth's mass inhomogeneities 
of third order are of smaller magnitude than previously supposed by most geodesists. 
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THE EQUATORIAL ELLIPTICITY OF THE EARTH 
FROM TWO MONTHS OF SYNCOM II DRIFT 
OVER THE CENTRAL PACIFIC 

by 

C. A. Wagner 

Goddard Space Flight Center 


INTRODUCTION 

This paper is one of a series dealing with the utilization of the Syncom II satellite to gain in- 
sight into the exact nature of the earth’s gravity field. During long periods between gross orbit 
change maneuvers initiated by on-board propulsion, this nearly synchronous satellite has remained 
over limited longitude regions of the earth for a time sufficient for it to accumulate noticeable 
drift perturbations caused by small longitude mass anomalies within the earth. The first gravity 
experiment with this satellite (Reference 1) was designed to test the sensitivity of its observed 
drift to the simplest kind of longitude- dependent gravity; namely, that associated with an earth 
whose equator is elliptical, not circular. As a result of that experiment the very strong assertion 
was made that the earth is, in fact, a triaxial ellipsoid whose gross elliptic equatorial ’’bulge” 
amounts to between 200 and 250 feet. This new gravity investigation, which essentially confirms 
the previous result, examines drift data from Syncom II for two months during the period April- 
June 1964. The mean daily drift rate of the satellite was 0.8 degrees/day during these months, and 
the nearly figure eight ground track progressively swept over 45° of longitude from 117° West to 
163° West. 

It was shown in Reference 1 that at a drift rate larger than 0.39 degrees/day a 24-hour satel- 
lite with the inclination of Syncom II would circulate completely around the world if left to drift 
freely. Clearly, Syncom II in April- June 1964 had a period far from "very close to synchronous.” 
The approximate theory of drift as a function of time for such a satellite, used in Reference 1 to 
analyse Syncom II data from August- February 1963-1964, was not serviceable for this "fast-drift” 
regime. The essential point is that the "very close to synchronous” theory dealt with long term drift 
in a limited longitude range. In fact, it assumed that the perturbation force averaged over the orbit 
referred to a single mean longitude position of the satellite, the longitude at synchronism. For 
data over a large longitude arc, it has been found most convenient and efficient to use a data re- 
duction theory stemming from the exact expression of drift rate as a function of longitude 
position for a 24-hour satellite. The important point in the applicability of this theory is that the 
net longitude drift of the satellite is limited over a single orbit. In the case of Syncom n during 
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April- June 1964, the small drift of 0.8° per orbit allowed the direct application of the theory to the 
ascending equator crossing of the satellite. For drift regimes with rates much larger than this, 
this simple theory for the drift of the ascending equator crossing probably should be modified to 
apply to the drift of the mean of the longitudes of the ascending equator crossings over a single 
orbit. Moreover, if the eccentricity of the 24-hour satellite is fairly large, R. R. Allan* suggests 
that this simple theory for the 24 -hour circular orbit satellite would still apply with the longitude 
position taken as the mean of ascending and descending equator crossings. 

At any rate, aside from the use of the first integral of the drift motion instead of the second, 
the analysis for longitude- dependent earth gravity in this report proceeds with the same basic as- 
sumptions as those of the previous report, namely: 

1. Only second-order longitude -dependent earth gravity (associated with an earth whose equa- 
tor is elliptical) was sensed during the long term drift of Syncom II in April-June 1964; 

2. All non-gravity perturbations on the orbit of Syncom n during this period were negligible; 

3. The earth zonal and sun and moon gravity effects on Syncom n drift can be adequately rep- 
resented over the span of a few months by a small bias in the reduced results. The bias 
can be estimated from simulated Syncom II trajectories, numerically integrated, as ex- 
plained below. 

Assumptions 2 and 3 proved good for the previous data analysis. Assumption 1 was not evaluated 
in that analysis. Calculations of probable higher order earth effects on the bias of recently re- 
ported geoids shows (Reference 2) that the Syncom n reduction in Reference 1 for the second-order 
effect only has a bias of no more than 25% in the reported magnitude of that effect. The major con- 
clusion of this report is that assumption 1 is probably satisfied to a greater degree than this for 
all the Syncom n drift data thus far available. 

The theory of Syncom II drift rate variation in a second order field is presented. The appli- 
cation of this theory to the calculation of the second order field from equator crossing data is then 
discussed. As shown in Reference 1, this second order drift theory is only strictly exact for a 
perfectly synchronous satellite, and ignores sun, moon and other earth perturbations. While the 
errors in the theory due to these sources were proved to be small in the previous Syncom analysis, 
there is no guarantee that they will be similarly small for the fast drift Syncom II in April-June 
1964. Therefore, in order both to verify the theory for this Syncom drift orbit and gain insight 
into the likely magnitude of these "model errors,” the field calculating method stemming from the 
drift rate theory is tested on simulated Syncom data. This simulated data consists of equator 
crossings closely paralleling the observed Syncom crossings for this period. 

The simulated data is derived by numerical integration of the motion of Syncom II starting 
from the observed elements in late April 1964. The closely paralleling numerically integrated 
Syncom trajectory is disturbed by the sun and moon. It is also disturbed by an earth which 


* Private communication. 
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includes, as inputs to the integrating program, two second order longitude gravity constants be- 
lieved to be close to those for the real earth. The proposed field calculating method, applied to 
this simulated data, estimates these input constants. The discrepancy between this estimate and 
the input constants is found to be small. This discrepancy is then taken as a good measure of the 
model error in the constants derived, in the final section of the report, from the observed Syncom 
crossing data. 


THEORY OF THE REDUCTION 

Equation 36 in Reference 1 gives the first integral of the drift motion of a 24- hour satellite in 
a second-order longitude- dependent earth gravity field as: 

(r ) 2 = A 22 cos 2y + Cj , ( 1 ) 


where y is the geographic longitude of the ascending equator crossing of the 24- hour satellite east 
of the equatorial minor axis, Cj is an integration constant for the motion, and A 22 is a function of 
the semimajor axis (a s ) and the inclination (i s ) of the satellite’s orbit, the mean radius (r o ) of the 
earth’s equator, and the magnitude (j 22 ) of second-order longitude gravity. A 22 is given by 


A 


22 



(cos 2 i s + l) 
2 


(2) 


Let A be the geographic longitude of the ascending equator crossing of the 24- hour satellite at a 
certain time. If /v 22 is the geographic longitude of the equatorial major axis, then K-y— ( 90 0 - A 22 ) 
from Figure 3 of Reference 2 and Equation 1 can be rewritten as: 


(A) 2 - ~A 22 cos 2 (A - \ 22 ) + Cj 


(3) 


It is evident from Equation 3 that if three sets of simultaneous values of A and A are available 
from position and velocity observations on the 24- hour satellite, then a unique set of values of A 22 , 
A 22 and C 1 can be calculated. These values, together with the known values of R 0 , a s , and i s , may 
be used in Equation 2 in order to calculate J 22 , the magnitude of the second-order longitude gravity 
effect. To show the form of the ’’condition equations” for A 22 and A 22 more explicitly, Equation 3 
may be expanded to give 


(A) 2 = Cj + cos 2A (-A 22 cos 2A 22 ) + sin 2A f-A 2 2 sin 2A 22 ) 


(4a) 
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or 


(A) 2 = Cj + c 2 cos 2A. + c 3 sin 2A 


(4b) 


where 


A. 


22 



(5) 


and 


A 


22 




1/2 


( 6 ) 


Thus, three or more sets of A and A values make it possible to determine a best set of C t , C 2 and 
C 3 values from Equation 4; for example, a best set in a "least squares" sense. 

From this set "best" values of A 22 and A 22 can be calculated from Equations 5 and 6. The 
actual observations on which the reduction from Equation 4 will be based are nine sets of orbital 
elements for Syncom II calculated at the Goddard Space Flight Center from range and range-rate 
radar data returned from the satellite during April- June 1964 (see Appendix A). For each of 
these elements, the longitude of the ascending equator crossing nearest to the epoch of the set plus 
six hours has been calculated numerically using a zonal gravity field through fourth order 
and a sun and moon gravity field, both of which are almost identical to the fields used in the orbit 
determination program for the satellite (Appendix A, Table A- 2). This particular ascending equa- 
tor crossing was chosen as that which fell as close as possible to the center of the observations 
from which the orbit in question was determined. From these equator crossing longitudes and 
times, longitude rates were estimated by the ratio of the change in longitude to the time elapsed 
between successive crossing points. In this way the relatively large daily- periodic effect of the 
sun and moon on the instantaneous "two-body orbit" drift rate was virtually eliminated from the 
"basic data" before the actual reduction for the second- order longitude gravity effect began. 

The question arises as to what geographic longitude, the longitude rate so calculated, belongs. 
If there were no acceleration in the longitude drift, it would make no difference where this calcu- 
lated drift rate was applied. But it is clear that with any small acceleration the correct 
application point is somewhere near the mean of the two successive longitudes. It is shown in 
Appendix C that approximating the exact longitude at which this calculated orbit rate should be ap- 
plied, by the mean of successive longitudes, is sufficiently accurate for the weekly data on Syncom 
II obtained in April- June 1964. 
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Table 1 


Drift Data Reduction from the Simulated Syncom II Trajectory of Appendix B. 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

i 

Time from 
Jan. 0.0, 1964 
(days) 

Mean 

Longitude, 

(degrees) 

Yi B (\) 2 

(“■* (S)1 

(degrees) 

X 2i = cos 2 

X 3i = sin 2^ 

Value of Y. 
Computed 
from C Given 
by Equation 15 

(”- (sgj) 

(»-• (^) ! ) 

e. 2 

(£?)') 

1 

120.6 

-120.49 

2.1005 

119.02 

-.48511 

.87445 

2.1003 

.02 

.0004 

2 

126.6 

-125.46 

2.0787 

109.08 

-.32689 

.94506 

2.0810 

- .23 

.0529 

3 

132.6 

-130.40 

2.0621 

99.20 

-.15988 

.98714 

2.0560 

.61 

.3721 

4 

138.6 

-135.31 

2.0254 

89.38 

.01082 

.99994 

2.0262 

- .08 

.0064 

5 

144.6 

-140.18 

1.9884 

79.64 

.17983 

.98370 

1.9928 

- .44 

.1936 

6 

150.6 

-145.01 

1.9566 

69.98 

.34235 

.93957 

1.9566 

.00 

.0000 

7 

156.6 

-149.78 

1.9090 

60.44 

.49302 

.86984 

1.9190 

-1.00 

1.0000 

8 

162.6 

-154.52 

1.8969 

50.96 

.63006 

.77671 

1.8807 

1.62 

2.6244 

9 

168.6 

-159.22 

1.8421 

41.56 

.74826 

.66340 

1.8430 

- .09 

.0081 

10 

174.6 

-163.86 

1.8029 

32.28 

.84545 

.53406 

1.8070 

- .41 

.1631 




^19.6626 


2.27791 

^8.57387 



4.4260 






2 

2^2.432133 

1 

2 

2^7.567804 






DATA REDUCTION OF THE SIMULATED SYNCOM II TRAJECTORY OF APRIL- 
JUNE 1964 FOR THE SECOND-ORDER LONGITUDE GRAVITY EFFECT 


This scheme of data reduction based on Equation 4 is now tested for a simulated Syncom II 
trajectory which begins with the actual orbit elements of epoch 64-04-25-02 U.T. in Appendix A 
and which closely approximates the actual Syncom II trajectory. In the digital program of this 
simulation (called "ITEM,” Interplanetary Trajectory by an Enke Method, at Goddard Space Flight 
Center), a basic earth zonal gravity and sun- moon field is used which is identical with that used 
in the actual Syncom II orbital determination program. In addition, a second-order longitude- 
dependent earth gravity field with J 22 - -1.68X10" 6 and\ 22 = -18° is used in the numerical 
computation of the two-month trajectory. (See Appendix B for the basic data of the simulated 
trajectory.) 

The drift data reduction for this simulated Syncom n trajectory— the reduction for the actual 
trajectory in the next section follows the same pattern— is described in what follows. Table 1 gives 
a summary of the reduction. 

Column 3 of Table 1 shows the mean longitudes, k iy upon which the data reduction is made. 
They are simply the mean longitudes of successive weekly equator crossings listed in column 3 of 
Table B-2. The corresponding drift rates, (column 4, Table B-2) are the differences in succes- 
sive longitude crossings divided by the time interval between these crossings. Column 4 of Table 
1 is self explanatory as are columns 5, 6 and 7. 


As stated in the previous section, the first problem is that of obtaining good values of c v C 2 
and C 3 (in Equation 4) from the data. The general procedure is as follows (this applies to the next 
section as well). Table 1 gives ten sets of corresponding values of (\) 2 and k, but they are incon- 
sistent in the sense that any three sets used in Equation 4 will give values of C lf C 2 and C 3 that in 
general are not the same as the values obtained for a different set of three. In order to make use 
of all the data, it is assumed that the k's are known exactly but that the values of the ( A .) 2 are sub- 
ject to slight random error. Then for the i th set (\.) 2 = c x +C 2 cos 2\. +c 3 sin 2k i + e. and the 
problem is that of establishing a criterion for estimating c l9 C 2 and c 3 from the ten equations i = 1 
to 10. The least squares method of estimation (Reference 3) seems most appropriate here, and it re- 
quires that those estimates of C v C 2 and c 3 be chosen which minimize the sum 



i = 1 


The simple change in notation indicated in the headings of columns 4, 6 and 7 of Table 1 is made 
in order to conform with the notation used in the statistical literature. When the change in notation 
is made, the equations take the form 


6 


(7) 


Y; 


= XnC, + X 2i C 2 


x 3 i c 3 + ^ 


10 , 


Least squares estimates of C 1? C 2 and C 3 (c l9 C 2 and C 3 ) are found and then the data in column 8, 
Table 1 is calculated from Y ci = X ti + X 2i C 2 + X 3i C 3 . Columns 9 and 10 of Table 1 are self- 
explanatory. 

It is shown in the theory of least squares estimation (Reference 3) that the c/s must satisfy 
the so-called normal equations: 



N 


£ v ‘ x “ 


N 


T Y i X 3i 






N 



1 


X 1 i X 2 



X 3i X 2i 



X 1 i X 3 i 




X 3 i X 3 i 


( 8 ) 


(9) 


( 10 ) 


From Table 1, those mixed- product sums for the simulated drift data which are not calculated at 
the bottoms of the columns in Table 1 are 


-4 /radV 

Y. X 2i — 4.06013 x 10 ydayy 


16.97526 x 10' 


^ x 2i x 3i = !• 


484560 . 


Equations 8, 9 and 10 for the simulated Syncom drift of Table 1 are therefore: 


19.6626 x 10' 4 = lO.OOOOCj + 2.277910, + 8.57387C 3 


(ID 


7 



4.06013 x 10“ 4 


16.97526 x 10 -4 



2. 2779 1C 1 + 2.432133C 2 + 1.484560C 3 


8. 57387CJ + 1.484560C 2 + 7.567804C 3 


(12) 


(13) 


The solution of Equations 11, 12 and 13 for C may be accomplished by any convenient algorithm. 
We choose to use the inverse of the right side coefficient matrix for this solution because the ele- 
ments of the inverse are also statistically valuable (see Equation 34). 

The inverse of the matrix of coefficients on the right sides of Equations 11, 12 and 13 is 
found to be 



8. 305427 

-2. 312087 

-8.955996 

(00EFF)' 1 = 

-2.312087 

1. 110737 

2.401565 


-8.955996 

2.401565 

9.807638 


When (COEFF) -1 is multiplied by the 1x3 matrix of left side constants of Equations 11, 12 and 13, 
the 1x3 matrix 



is found to be 


C(simulated data) 




/ rad\ 2 

1.88855 x 

10~ 4 

(agj 

/ rad\ 2 

-0.18471 x 

10' 4 

[day) 



( rad\ 2 

0. 13970 x 

10“ 4 

yd ay J 


(15) 


From Equations 15 and 5 the "best" value of ^ 22 for the drift data in Table 1 is 


k - o’ tan 

22 , s im 1 


-c, 

-c. 


-0.13970 


2 tan " 1 0.18471 


*18.4° 


(15a) 
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Similarly, Equations 15 and 6 give 


(^ 22 ) simulated = s im . + C 3 % im .) 1/2 = 23 . 15 X 10' 6 ' 


23. 15 x 10" 6 x 


= 23.02 x 10~ 6 


/l436. 06817 solar day\ : 
\1440. 00000 sid. day/ 

/ rad y 
^sid. day/ 


(16) 


Data from Table B-l, Appendix B used in Equation 2 yields 


, v V r o / 6378.4 V 

( A 22 ) s i m . “ (J 2 2)s im . 72 x 77 x 1^42230. 01 J 

( J 22 ) Sim. (“I 3 - 86) (sid. day) 


(cos 2 32. 6° + 




(17) 


From Equations 16 and 17 the "best" value of J 22 for this simulation, according to the drift theory 
presented in Chapter 2, is 


J 


22.sim. 


23.02 x 10~ 6 
13.86 


-1.66 x 10' 6 


(18) 


The simulation actually used the following values of J 22 and \ 22 : 

^ 22 < sim . = -1-68 X 10- 6 , 

^"22. sim. = -18 0° • < 19 ) 

Comparison of Equations 15a and 18 with Equation 19 shows that the ,T model error” of the drift 
theory presented in the previous section, which ignored drift rate changes due to sun, moon and 
earth zonal gravity influences, accumulates a small bias in the results for J 22 and^ 22 calculated 
according to the simple theory for a two- month simulated Syncom II drift which closely approxi- 
mates the actual drift in this period (Figure 1). This model bias (ignoring higher order earth 
longitude gravity effects) in J 22 calculated from the theory for this two- month drift is estimated 
from the above results of the simulated data reduction as 

J 22 ,s im . - J 22 .sin, = -1-68 X 10- « + 1.66 - 10“ 6 


= -0.02 x 10‘ 6 


( 20 ) 


9 



DRIFT RATE: - °/ DAY 


18.0 + 18.4 


+ 0 . 4 ' 


( 21 ) 


The model bias in \ 22 is estimated as 

^22»sim. ^22.sijn. 


REDUCTION OF NINE SYNCOM II ORBITS IN APRIL-JUNE 1964 TO DETERMINE THE 
ELUPTICITY OF THE EARTH'S EQUATOR 


Table A-l in Appendix A lists nine orbits in April- June 1964 for Syncom II computed at GSFC 
from range and range- rate information without taking into account the anomalous gravity field due 
to the ellipticity of the earth's equator. Figure 1 shows the observed steady decrease in the west- 
ward drift rate of the mean longitude of the satellite for 45° of drift excursion over the Central Pa- 


cific in this period. This rate decrease is closely 



MEAN LONGITUDE (DEGREES WEST) 


paralleled by that in the simulated drift of Syncom 
II using a longitude- dependent gravity field 
which includes sun, moon and high- order earth 
zonal components. Other simulated trajectories 
have shown minor, negligible drift rate changes 
over a duration of time measured by months if 
longitude -dependent earth gravity is not used in 
the simulation. The experience with the simu- 
lated Syncom II trajectory (the previous section 
and Appendix B), which closely agrees with or 
closely parallels the actual drift, gives us con- 
fidence that the simple reduction theory used 
successfully in the previous section to deter- 
mine equatorial ellipticity from the simulated 


A. NUMERICALLY CALCULATED DRIFT RATE 
IN A PARTICLE TRAJECTORY FROM THE 
SYNCOM 2 ELEMENTS OF EPOCH: 64 - 04 - 25 - 2.0 HRS, 
WITHOUT LONGITUDE-DEPENDENT EARTH GRAVITY 

o DATA FROM A SIMULATED SYNCOM U TRAJECTORY 
BEGINNING WITH THE ELEMENTS OF EPOCH 
64-04-25-02; USING A SUN, MOON, ZONAL EARTH 
& LONGITUDE-DEPENDENT EARTH GRAVITY FIELD 
WITH: 

J 22 - -1 .68X10” 6 
X 22 - 18° 


data can be applied as successfully to determine 
equatorial ellipticity from the actual data. The 
following Table 2 of drift data has been derived 
from Appendix A using the nine orbits of Syn- 
com II in April-June 1964. 

From Table 2, the mixed product sums not 
calculated at the bottoms of the columns are: 


“ THEORETICAL DRIFT IN A SECOND ORDER LONGITUDE- 
DEPENDENT EARTH GRAVITY FIELD ONLY WITH: 

J 22 - - 1 .66X10" 6 
A 22 =-18.4° 

• DATA FROM ACTUAL SYNCOM II OBSERVATIONS 

— THEORETICAL DRIFT IN A SECOND ORDER LONGITUDE- 
DEPENDENT EARTH GRAVITY FIELD WITH: 

J 22 = -1 .69X10 “ 6 
\22 = -17.5° 


NOTE: THE DISPLACEMENT OF THE TWO SETS OF DATA 

IS DUE MAINLY TO THE ERROR IN THE SEMI'MAJOR 
AXIS FOR SYNCOM II, REPORTED FOR EPOCH 
64-04-25-02 



2.87010 x 10“ 4 




13. 17124 x 10" 4 



x 2 i x 3i 


Figure 1— Drift rate as a function of mean longitude for 
Syncom II during April-July 1964. 
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Table 2 


Drift Data Reduction for Syncom II, April- June 1964. 


Mean . 1 

Longitude (\.) 2 2^ X 2i 

K ( /rad\ 2 \ (degrees) cos 2\ £ 

(degrees) ^10" 4 ^dayj j 


X 


3 i 


sin 2Xi 


Y C i 

(Using Cfrom 
Equation 26) 




1 

-121,22* 

2.004 

117.55 

-.46252 

.88661 

2.002 

.20 

.040 

2 

-128.09 

1.973 

103.81 

-.23870 

.97109 

1.970 

.30 

.090 

3 

-133.70 

1.934 

92.61 

-.04553 

.99896 

1.936 

- .20 

.040 

4 

-138.84 

1.882 

82.32 

.13364 

.99103 

1.9005 

-1.85 

3.423 

5 

-144.72 

1.868 

70.56 

.33282 

.94299 

1.8535 

1.45 

2.103 

6 

-150.935 

1.810 

58.13 

.52799 

.84925 

1.805 

.50 

.250 

7 

-155.94 

1.770 

48.13 

.66744 

.74466 

1.760 

1.00 

1.000 

8 

-160.83 1 ' 

1.716 

38.35 

.78424 

.62046 

1.720 

- .40 

.160 



14.957 


^ 1.69938 

2 ^ 7.00505 



22 7.106 





2^ 1.74089 

2 

2^6.25910 





♦On 30 April 1964. 
tOn 20 June 1964. 



Following the technique of data reduction according to the 24-hour satellite drift theory dis- 
cussed in the previous section, the normal equations for the "best" values (in the least squares 
sense) of the ellipticity-drift rate parameters c are: 


4 /rad\ 2 
14.957 x 10- 4 

= 8Cj + 1. 69938C 2 + 7.00505C 3 , 

(22) 

4 /rad \ 2 

2.87010 x 10" 4 ^day/ = 

1. 69938Cj + 1. 74089C 2 + 1.19093C 3 , 

(23) 

. /rad\ 2 

13. 17124 x 10" 4 (day^J = 

7. 00505Cj + 1.19093C 2 + 6.25910C 3 . 

(24) 


The inverse of the matrix of coefficients on the right sides of Equations 22, 23 and 24 is: 



14.008120 

-3.3905080 

-15.032469 

(COEFF)" 1 = 

- 3.3905080 

1.4810135 

3.5127885 


-15.032469 

3.5127885 

16.315400 


The ellipticity-drift rate parameters 


(Q 

(Q 

(C 3 ) 


are then found to be: 


C( actual data-unadjusted for model error) 


fr ad\ 2 

1.79210 x 10~ 4 (di^J 

f ra d V 

-0.19339 x 10- 4 

4 ( ra d \ 2 
0.13546 - 10- 4 [^) 


(25) 


(26) 


Thus, from Equations 26 and 5 the best value (unadjusted) of a. 22 (actual) as sensed by the drift of 
Syncom n in April-July 1964 is 


X- 22 (actual data-unadjusted) 



1 

2 


tan" 1 


/-0. 135 46\ 
l 0. 19339 J 


= -17.5° 


(27) 
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For the Syneom II elements in Appendix A, using a s and i s from the second set of elements in the 
drift period (since they are probably better determined than the first), Equation 2 gives 


A 


22 


J22 


-72t 7 2 x 


^6378.4 \ 2 ~[ ["cos 2 32.6 ° + 1 
42228. 8 J |_ 2 


= "J 22 (13.859) 



(28) 


From Equations 26 and 6, the best value (unadjusted) of A 22 (actual data) as sensed by the drift of 
Syneom II in April-July 1964 is 


A 22 ( actu al data -unadjusted) 


(c/ + c/) 


1/2 


23.6 x 10' 6 


r ad 


Isolar 



2 


= 23.6 x 10“ 6 


( r ad A 2 

^solar day / 


x 0.99455 


/ sol . 
^sid. 


day\ 2 
day/ 


= 23.47 x 10“ 6 


( 


rad \ 
sid day/ 


2 


(29) 


Equations 28 and 29 show that the best value (unadjusted) of the magnitude of the equatorial 
ellipticity of the earth sensed by Syneom II in April-July 1964 is 

J 22 ( actual dat a-unadj usted) ” -1.69 x 10 -6 (30) 

In adjusting the values given by Equations 27 and 30 for T, model error" by means of the biases 
found in the simulated Syneom II trajectory, it is assumed that all influences on Syneom II drift in 
this period other than earth, sun and moon gravity, are negligible (see DISCUSSION). The data 
reduction of the simulated Syneom II trajectory in April- June 1964 in the previous section showed 
that the simple drift theory presented herein produced a bias in the reduced J 22 of -.02 x 10 ~ 6 and 
a bias in the reduced \ 22 of +0.4°. Application of these biases adjusts Equations 27 and 30 for all 
but higher order longitude -dependent earth gravity effects and gives 

J 22 ( actual data-adjusted) = -(1.69) x 10' 6 “0.02 x 10" 6 = -1.71 x 10~ 6 , ( 31 ) 

\ 22 (actual data-adjusted) - “17.5° +0.4° = “17.1° . (32) 

To gain insight into the allowable variations of J 22 and ^ 22 due to the errors of observation 
and effects not accounted for (model errors) in this gravity experiment, the simplifying assumption 
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II 


is made that all deviations in Tables 1 and 2) of the actual or simulated data from the theoreti- 
cal data are random and are normally distributed with mean zero and variance cr 2 . It is to be 
noted that the "best values" from the simulated data reduction established a most likely "system- 
atic" error in J 22 and \ 22 for the actual experiment. This systematic error was due to all im- 
perfections in the simple model tested by the simulation and equivalent imperfections are likely to 
be present in the real data (see DISCUSSION). 

If the above assumption is made for the deviations e . in the simulated data, then an unbiased 
estimate of cr 2 (c sim ) calculated from the data of Table 1 is 


S 2 


4.4260 x 10" 12 / rad\ 4 
10-3 \da y y 


0.6323 x 10' 12 



An unbiased estimate of the standard deviation S is 


S(« sim .) = (0.6323x 10* 12 ) 1/2 = 0.7952 x 10- 6 • (33) 

It can be seen that only for i = 8 in Table 1 is the error greater than 2S(e sim J for the experiment 
as a whole. 

The statistical theory for normal distributions (Reference 3) shows that the parameters C } of 
the normal equations have expected values C. (the actual drift- ellipticity parameters tested for) 
and variances a} , where 


cr 2 C j = (cr 2 ( £ )) (COEFF j j ) ' 1 • 

(oOEFFj j ) “ 1 denotes the diagonal element in the ith row and jth column of the inverse of the coef- 
ficient matrix of the normal equations for the experiment. An unbiased estimate of is 

C: 


S 2 = (s 2 ( £) ) (COEFF..)-! . 

It follows from Equation 34, using the values given by Equations 33 and 14, that 


(34) 


Sp = 0.7952 Vl. Ill' x 10' 6 

2, s im. 


Sg = 0.7952 V9.808' x 10' 6 

3, sira. 


0.0084 x 10- 4 (^) 2 ’ 
a /rad\ 2 

0.0250 x 10' 4 (a^) 


(35) 


(36) 
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Thus the standard deviation bounds on C 2 sim and C 3 sim are 

C 2 ,s im . = -(0.1847 ±0.0084) X 10- 4 (g) 2 - 

^3, sim. = (0.140 +0.025) x 10' 4 

The same theory applied to the actual Syncom II data in Table 2 yields 

x 7. 106 x 10~ 12 / rad\ 4 10 /rad \ 4 

S2 ( e act.) = ~8^3 id^j = 1.4212 x 10' 12 ^dayj 

and 

s ( e act.) = 1. 192 x 10“ 6 (iJay) 2 * (39) 

Comparison of Equations 39 and 33 shows that the standard error of the simulated experiment 
(which is free from "observation error") is not much smaller than that of the experiment on the 
actual data. In Reference 1 it was found that the actual experiment had a standard error about five 
times that of the simulated. Evidently the set of orbits for Syncom II in April- June 1964 is of 
much greater precision than those of the two earlier long drift periods. 

For the actual Syncom II data treated in the same fashion as the simulated data, the gravity 
experiment leads to the bounds on C computed below. 

Unbiased estimates of ag and og are 


(37) 

(38) 


and 



1. 192 /nr x io~ 6 


= 0.0145 x 



Sp = 1.192 /16.3154' x 10~ 6 = 0.0481 x 10" 4 

^3, act. 

Thus the standard bounds on C 2 and C 3 in the actual experiment are 

l rad\ 2 

S 6 - -(0.1934 ±0.0145) x 10" 4 JiW 

2 , actual data \ y / 

. /rad\ 2 

Sp = (0.1355 +0.0481)x 10~ 4 dav 

3, actual data \ y * 



(40) 

(41) 


(42) 

(43) 
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It can be seen that though the standard error of the actual experiment is only 1.5 times that of 
the simulated, the standard bounds on C are more nearly twice that in the simulated experiment. 
This is most probably due to the fact that two additional data points which sense the equatorial el- 
lipticity were used in the simulated experiment. 

In establishing standard bounds on f 22 and k 22 from the experiments, all possible combinations 
of these quantities from the standard bounds on C 2 and C 3 are not computed. Clearly, the chance of 
obtaining c 2 + and c 3 + simultaneously from the experiment is less than the chance of ob- 
taining a single standard deviation for each separately. In Appendix E the bounds on the adjusted 
values of J 22 and ^ 22 sensed by the Syncom II drift in April- July 1964 are computed as 


J 22 (actual data) = -( 1 . 71 ± 0-22) x 10“ 6 


(44) 


and 


k 2 2 (actual data) - -(17.1 ±4.9°) 


(45) 


DISCUSSION 

The Syncom II drift data reduction for sensitivity to the equatorial ellipticity of the earth, used 
in this report, is open to a number of lines of criticism. 

In the first place it may be fairly asked whether the observed decrease in drift rate can be ex- 
plained by some force not considered here. \The direct evidence of agreement in almost all partic- 
ulars between the observed trajectory of Syncom II during the period analysed and a simulated 
trajectory using a recent set of gravity constants for the earth, sun and moon make it unlikely that 
such an '’outside force” could have had appreciable influence on the drift. A number of such pos- 
sible outside influences on the orbital energy of the satellite were discussed in Reference 1. Solar 
radiation pressure was found to have an effect on Syncom II less by many orders of magnitude than 
that of solar gravity. It has never been proven that gas leaks from the on-board propulsion units 
during the long drift periods. Even if it did, it would not change the energy of the satellite appre- 
ciably provided that the accidental unrecorded ”on times” of this leakage were reasonably distributed 
around the spin cycle of the satellite. The magnitude of effects due to micro meteorite drag on the 
satellite is difficult to estimate because of insufficient data on this phenomenon. 

With respect to all "outside influences” except higher order earth longitude gravity, the almost 
exact prediction on the basis of longitude- dependent earth gravity alone of the long term drift de- 
celeration between 120° West and 165° West from the observed acceleration between 55° West and 
65° West (Reference 1) makes it even more unlikely that any such influence had appreciable effect 
on Syncom n during its "free drift” periods. 
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A second line of criticism of the gravity results summarized in Equations 44 and 45 stems 
from the uncertainty of the separation of systematic and random errors in the experiment. The 
assumption made in the statement of error bounds in Equations 44 and 45 is that all errors in the 
observed drift rate are randomly distributed about that theoretically calculated from the model. 

If the unknown "systematic error” due to high-order longitude gravity is ignored, the results of the 
simulated trajectory reduction appear to show that "model error" was the principal "noise" source 
in the actual gravity experiment. This is to say that the systematic error due to the insufficiency 
of the model to predict the exact results in a full gravity field gives, when treated as a random 
effect, a standard experimental error only somewhat smaller than that associated with the actual 
data. Furthermore, one can see directly from Figure 1 the close parallelism of drift rate devia- 
tions between the actual and simulated trajectories. These facts suggest that the small overall 
systematic errors of the simulated experiment (for example, .02 x 10 -6 in J 22 and 0.4° in k 22 ) are 
closer to the actual order of magnitude of the standard errors of the ellipticity parameters than 
that reported in Equations 44 and 45. In other words, the "best" results from the actual experi- 
ments are perhaps closer to reality (neglecting higher order gravity effects), with bounds half 
as wide as reported and with the same degree of confidence. 

On the other hand, preliminary, less refined reductions of more limited Syncom II data in this 
period gave evidence of a longitude effect near the limits of the bounds in Equations 44 and 45 with 
only a small increase in the standard experimental error. Hence it is felt that until a more refined 
analysis of this data is undertaken to remove "model error" point for point, Equations 44 and 45 
give fair, at best conservative bounds on the ellipticity of the earth's equator as seen by Syncom 
II during two months of drift over the Central Pacific. 

When these results are compared with those for Syncom II drift over Brazil (Reference 1) and 
with expected results on the basis of higher order effects as predicted from other sources of grav- 
ity data (Reference 2), it would appear that when the full longitude field is finally mapped the "true" 
value of J 22 will be closer to -1.7 x 10~ 6 than -1.9 x 10“ 6 as predicted in Reference 2. What this 
means is that the earth's third- order longitude mass deviations (to which Syncom II is next most 
sensitive) are probably of smaller extent than previously supposed. Exact statements on the bounds 
of higher order gravity effects as sensed by Syncom II must wait on an integrated analysis of all the 
drift data for this satellite. 


CONCLUSIONS 

1. The steadily decreasing drift rate of Syncom II during a two- month excursion over the Cen- 
tral Pacific in April- July 1964 confirms the extent and location of earth equatorial ellipticity pre- 
dicted earlier from observations on this satellite over Brazil. 

2. The results of this new Syncom n gravity experiment show a difference in the major and 
minor radii of the earth's elliptical equator of 

a Q -b 0 = 214 ± 28 feet jj 22 o -(1. 71 ±0.22) x 10” 6 ] . 
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The longitude location of the semi- major axis of the earth’s elliptical equator as sensed by Syncom 
II during this period is 


k 22 - -(17.1 ±4.9)° . 

3. The above results combined with the results of the independent gravity experiment with this 
satellite when it was over Brazil (Reference 1) indicate that the earth is more homogeneous, to at 
least third order, than has been supposed. 

4. Tentatively, the M best parameters” of equatorial elliptic ity as seen by Syncom n to date are 

J 22 = -1.7 x 10~ 6 , 

\ 22 = -18° . 

5. Only second- order longitude- dependent earth gravity need be considered in the design of 
station keeping systems for future 24 -hour satellites, and this without appreciable modification for 
longitude location. 
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Appendix A 


Orbital Elements and Partially Reduced Drift Data for 
Syncom II During April-June 1964 


Table A-l gives the elements for Syncom n during April- June 1964 as derived from radar 
range and range-rate data in the Goddard Space Flight Center T s Data Systems and Tracking Direc- 
torate. By means of the particle trajectory digital program called "Item" at GSFC, ascending 
equator crossing longitudes and times nearest to the epoch dates in Table A-l plus six hours were 
calculated. These are listed in Table A-2. The choice of epoch date plus six hours was mainly 
dictated by the desire to use the crossing that was closest to the "center of gravity" of the data 
that went into the determination of the orbit. The gravity field used for the short term trajectories 
in these calculations was virtually the same as that used in the orbit determination program that 
calculated the elements used at the start of each trajectory (see Appendix D). The only difference 
was the inclusion of J 22 and \ 22 in the earth potential of the particle trajectories. Since the values 
for these constants used in the "Item" program are close to "reality," it is likely that the ascend- 
ing equator crossings so calculated are even more accurate than an estimation of them made by 
using precisely the same field as that in the orbit determination program. 

Table A-l 


Orbital Elements for Syncom II in April-June 1964 as Calculated in the Data Systems and Tracking 
Directorate at the Goddard Space Flight Center from Radar Range and Range Rate Data (Orbit 12). 


Crossing 

Number, 

j 

Epoch (yr-mo- 
day-hr-min, U.T.) 

Semimajor 
Axis, a s 
(km) 

Eccentricity, 

e 

Inclination, 

i 

s 

(degrees) 

Argument of 
Perigee, o» 
(degrees) 

Mean 

Anomaly, 

m 

(degrees) 

Right 
Ascension 
of the 
Ascending 
Node, B 
(degrees) 

1 

64-04-25-02-00.00 

42230.01 

.00119 

32.603 

198.716 

333.752 

313.879 

2 

64-05-05-16-00.00 

42228.81 

.00123 

32.562 

199.105 

185.331 

313.739 

3 

64-05-12-16-00.00 

42228.04 

.00127 

32.625 

201.959 

183.885 

313.592 

4 

64-05-19-14-00.00 

42227.80 

.00119 

32.577 

199.859 

157.342 

313.548 

5 

64-05-25-15-00.00 

42227.65 

.00128 

32.599 

200.270 

173.279 

313.401 

6 

64-06-02-21-00.00 

42226.72 

.00130 

32.578 

201.136 

264.160 

313.331 

7 

64-06-09-21-00.00 

42225.74 

.00123 

32.579 

199.640 

267.215 

313.247 

8 

64-06-16-15-00.00 

42224.41 

.00122 

32.567 

200.824 

177.655 

313.141 

9 

64-06-23-15-00.00 

42224.17 

.00129 

32.561 

200.344 

179.909 

313.020 
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Table A-2 


Drift Data for Syncom II During April- June 1964 (Orbit 12). 


1 

2 

3 

4 

5 



Longitude of the 



Crossing 

Number, 

j 

Time from 
Jan. 0.0, 1964 
(days) 

Ascending Equator 
Crossing Nearest 
to Orbit Epoch 
+6 Hours 

K 

(degree s/day) 

d> 2 

(”-* i^y) 



(degrees) 



1 

116.6015 

-117.169 

.8111 

2.004 

2 

126.5986 

-125.278 

.8047 

1.973 

3 i 

133.5948 

-130.908 

.7969 

1.934 

4 

140.5910 

-136.483 

.7860 

1.882 

5 

146.5871 

-141.196 

.7830 

1.868 

6 

155.5820 

-148.239 

.7708 

1.810 

7 

162.5776 

-153.631 

.7622 

1.770 

8 

168.5735 

-158.201 

.7505 

1.716 

9 

175.5687 

-163.451 




NOTES: 


The data in columns 2 and 3 were derived from the elements in Table A-l by the numerical calculation of particle trajectories 
in an earth, sun and moon gravity field specified by the constants in Table B-2. The longitude rates in column 4 are the first 
longitude differences of column 3 divided by the first time differences of column 2. 
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Appendix B 


Orbital Elements and Partially Reduced Drift Data from the Simulated 
Syncom II Trajectory During April-June 1964 


Table B-l gives the set of elements at weekly ascending equator crossing times during a sim- 
ulated Syncom II trajectory (computed numerically by GSFC T s "Item" Program) beginning with 
the actual elements of Syncom II, epoch 64-04-25-02. The gravity field of the particle trajectory 
was exactly that used in the orbit determination program for Syncom II with the addition of equa- 
torial ellipticity introduced into the geoid of amount and location given by the constants 

J 22 = “1.68 x 10“ 6 


and 


X 22 = ”18° - 

(See Appendix D for a representation of the geoid used in the simulation.) 


Table B-l 

Orbital Elements For the Simulated Syncom II Trajectory Beginning with 
the Syncom II Elements of Epoch 64-04-25-02-00.00 UT (Orbit 12). 


Crossing 

Number, 

j 

Semimajor 
Axis, a s 
(km) 

Eccentricity, 

e 

Inclination, 

i 

s 

(degrees) 

Argument of 
Perigee, 
(degrees) 

Mean 

Anomaly, 

m 

(degrees) 

Right 
Ascension 
of the 
Ascending 
Node 
(degrees) 

Epoch 
(days from 
Jan. 0.0, 1964) 

1 

42229.86 

.00121 

32.601 

196.567 

163.665 

313.851 

117.604 

2 

42230.60 

.00123 

32.601 

199.401 

161.849 

313.758 

123.601 

3 

42229.94 

.00123 

32.585 

197.459 

164.838 

313.675 

129.598 

4 

42229.31 

.00122 

32.583 

200.538 

162.840 

313.567 

135.595 

5 

42229.30 

.00116 

32.570 

197.176 

163.552 

313.490 

141.592 

6 

42227.75 

.00122 

32.566 

198.568 

163.326 

313.388 

147.589 

7 

42228.87 

.00119 

32.562 

200.113 

162.958 

313.317 

153.585 

8 

42226.32 

.00126 

32.549 

198.869 

161.706 

313.220 

159.582 

9 

42227.77 

.00113 

32.547 

200.058 

161.779 

313.142 

165.578 

10 

42225.77 

.00117 

32.532 

197.797 

165.382 

313.055 

171.575 

11 

42226.02 

.00120 

32.533 

200.697 

160.098 

312.972 

177.571 
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In Figure B-l the elements numerically calculated during this trajectory (semi- major axis and 
inclination) are compared with those actually reported for Syncomll in April- June 1964. It is noted 
that among all the elements, only the semimajor axis shows a significant discrepancy in the simu- 
lated trajectory. The results in the two sections on data reduction show that the long term trend of 
the semi- major axis changes after the initial epoch is very well explained by the drift theory 
of this report. It should be understood that instantaneous drift rate changes of the ground 
track of the orbit of the 24-hour satellite are coupled directly to semimajor axis changes by 
Kepler’s third law which relates the semimajor axis to the period of the instantaneous orbit and 
the central, or gaussian, gravity constant of the principal attracting body. 

It is seen in Figure B-l that the simulated semimajor axes are on an average about one kil- 
ometer greater than the actual ones throughout the trajectory. This agrees well (in terms of 
Kepler’s third law) with the drift rate differences between the two trajectories at corresponding 



THE SIMULATION USED AN EARTH FIELD WITH THE 
SECOND ORDER LONGITUDE GRAVITY CONSTANTS: 
J 22 = -1-68X KT 6 
*22 =- 18 ° 

THE DOWNWARD TRENDS OF THE SEMIMAJOR AXIS 
DATA IN BOTH ACTUAL AND SIMULATED CASES ARE 
CLOSELY PARALLEL. 


THE APPARENT "DISPLACEMENT" IN THE TWO SETS OF DATA 
IS PROBABLY DUE TO ERROR IN THE SEMIMAJOR AXIS OF 
SYNCOM H AT EPOCH 64-04-25-02 U.T. 

NOTE: • ACTUAL SYNCOM II DATA 

O SIMULATED SYNCOM H DATA ( FROM THE 
NUMERICALLY INTEGRATED TRAJECTORY) 


Figure B-l— Orbital data from S/ncom li observations in April-July 1964 compared to data from a numerically 
integrated trajectory in a ful ! earth , sun and moon field, beginning with the Syncom II elements of Epoch 64-04-25-02 U.T. 



times (Figure 1). This discrepancy then is probably due to the relatively large error in the deter- 
mination of the semimajor axis for Syncom II (and used in the simulation) at the start of the drilt 
period, epoch 64-04-25-02. 

Table B-2 gives drift rate and longitude data from the numerically computed simulated trajectory. 


Table B-2 


1 

Crossing 

Number, 

j 

1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 


Data From a Two -Month Simulated Syncom II Trajectory Beginning with 
the Syncom II Elements of Epoch 64-04-25-02-00.00 UT (orbit 12). 


2 


Time from 
Jan. 0.0, 1964 
(days) 

117.6034 

123.6006 

129.5977 

135.5946 

141.5916 

147.5884 

153.5851 

159.5816 

165.5781 

171.5744 

177.5706 


3 

Longitude at 
the Ascending 
Equator Crossing 
(degrees) 

-118.000 
-122.980 
-127.934 
-132.868 
-137.758 
-142.603 
-147 .409 
-152.156 
-156.888 
-161.551 
-166.164 


4 

k 

(degrees/day) 

.83039 

.82607 

.82276 

.81541 

.80793 

.80144 

.79163 

.78913 

.77765 

.76932 


5 



2.1005 

2.0787 

2.0621 

2.0254 

1.9884 

1.9566 

1.9090 

1.8969 

1.8421 

1.8029 


NOTES: Initial Elements and Gravity Constants of the Trajectory: 


semimajor axis, a = 42230.01 km 
inclination, i = 32.603 degrees 
mean anomaly, m = 333.752 degrees 
mean radius of equator, R 0 = 6378.388 km 
gravity constant, M sun = 332,490 dearth 
angular velocity, co e = 0.72921 1 5 x 10~ 4 rad/sec 


20 

= 1 . 08219 x 10 3 

J 30 

-2 . 285 

40 

- -2.12 x 10" 6 

^22 

-1.68 x 


eccentricity, e =0-00119 

argument of perigee, = 198.716 degrees 

right ascension of the ascending node = 313-879 degrees 

gravity constant, dearth = 3-98627 x 10 5 km 3 /sec 2 

gravity constant, M moo n -0. 01229491 Mearth 

x 10“ 6 

i°“ 6 ^ 22 = " 18 . 0 ° - 
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Appendix C 

The Error in Estimating that the Mean Velocity of an Accelerating Trajectory 
Between Two Points in Space Occurs Midway Between the Points 

Suppose that a point m at t = 0 is at S = S 0 with velocity v - v 0 and moves with constant 
acceleration, a. On the assumption that S, v and a are colinear, the (s, t) graph for the motion 
is given by 

S = | at 2 + v 0 t + s 0 . (Cl) 

The graph is plotted in Figure Cl. The velocity, v, at any time in the interval is 


ds 

dt 


at + v 0 


(C2) 


and the mean or average velocity, v, over the time interval 1 1 is 


v o + V tl 


(C3) 


since the motion is uniformly accelerated. But 
at 1 1 the velocity, from Equation C2 is 


at i + v o 


(C4) 


Hence, from Equations C3 and C4, 


(C5) 


Comparison of Equations C5 and C2 shows that 
v is the instantaneous velocity of the point at the 
t = t j/2 . The position of the point at this time 
is, from Equation Cl, 



S' 


1 AiY /V\ 

2 a \2j + v o Y 2 ) 


+ S n 


Figure C-l . The (s, t) graph corresponding to the one- 
(C6) dimensional trajectory of a point with constant acceleration. 
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Now As (Figure C-l) is clearly 


As " 2 at i 2 + v o * 1 = v A t . 


(C7) 


Hence in estimating s'by s Q +As/2 an error As e is introduced given by 


As z i 

As £ - s 0 + 2 “^ ” S Q + a 4+v 


o 2 


8 ' v o 2 


o “ S. 


at i _ at l 

4 8 


or 


As c 


at 1 
8 


(C8) 


If a or t x 2 or both are small, As may be sufficiently small to validate the estimation. 


Syncom H drift is given, in terms of the longitude acceleration, as 


A- - A 22 sin 2 (K~ ^22) (^ 9 ) 

Considering the right side of Equation C9 a constant for a limited excursion of X, the solution of 
Equation C 9 is in the form of Equation Cl with 

K - s 


and 


A 22 sin2(^-\ 22 ) - a. 


The maximum value of |A 22 sin 2 (a. - k 22 ) | is reached when 2 (a. - A.) = ±77/2 , and +377/2. From 
the earlier Syncom II gravity experiment (Reference 1), it was determined that J 22 = -1.7 x 10“ 6 . 
For a s = 42230 km and i s = 32.6° (the Syncom n orbit parameters in April- June 1964), 


A 


2 2 


(cos 2 i 

-72" 2 (R 0 /a s ) 2 y 


1 ) 


-23.0 x 



Thus | a (max. )| for this drift period is 


/rad \ 2 

23.0 x 10- • • 



From Equation C8, the maximum error made in assigning the mean interval drift rate of Syncom H 
to the mid- longitude of an interval traversed in ten days is 


AX 


e 


23 x 1(T 6 x 100 x 57.3 
8 


0.0165° . 


This error is sufficiently small to be ignored in the overall data reduction. Thus, in assigning the 
mean longitude in a Syncom II drift interval of the order of 10 days to the mean drift rate in that 
interval, one need retain only two decimal figures in the longitude location. 
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Appendix D 


The Earth Gravity Potential and Force Field Used in this Report 

The gravity potential used as the basis for the data reduction in this study is the ex- 
terior potential of the earth derived in Reference D-l for geocentric spherical coordinates 
referenced to the spin axis and the center of mass of the earth. The infinite series of 
spherical harmonics is truncated after J 44 . The zonal harmonic constants J 20 , J 30 and J 40 
are derived from Reference D-2. 

The radius of the earth, R 0 , used in this study is 

R 0 11 6378.388 km 

The earth’s gaussian gravity constant used is 

km 3 

= 3.9862677 x 10 s ~ * 

a sec 

Neither of these values, taken from Reference D-3, nor the ’’zonal geoid” of Reference D-2, is felt 
to be the most accurate known to date. They are the values used by the GSFC Tracking and Data 
Systems Directorate to calculate the orbit elements of Syncom II from radar and Minitrack obser- 
vations. They were chosen to insure consistency between the data of this study and these published 
orbits, inasmuch as the "triaxial” reduction for which this study has been undertaken is not signifi- 
cantly affected by the probable errors in these values. The second-order longitude harmonic con- 
stants used in the simulation studies were 

J 22 - -1.68 x 10“ 6 

and 

These are the values shown or implied on the ’’longitude geoid’’ below (for the J 22 harmonic). At a later 
point in the analysis, the slightly different values reported in the abstract were estimated. The most 
accurate ’’zonal geoid” is probably that of Kozai (cited in Reference D-4) with the following 
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earth constants 


R 0 = 6378.2 km , J 20 = 1082.48 x 10“ 6 , J 30 = “2.56 x 10“ 6 


= 3.98603 x 10 s 


km 3 


-1.84 x 10" 


The earth T s gravity potential as used in the simulations of this report (to fourth order, probably 
sufficient to account for all significant perturbations on a 24- hour satellite) may be illustrated 
(following Reference D-4, Appendix B, with the zonal constants of Reference D-2) as follows: 


NORTH 




J30 Rq 3 

2r 3 


(5 sin 3 <j> - 3 sin <£ ) 


P 3 
K o 


2r 3 


cos 4> (l5 sin 2 4- - 3 ) cos (k - \ 31 ) 
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8r 4 (35 sin 30 sin 2 <£+ 3 ) - (l40 sin 3 0 - 60 sin tf>) cos cos (\ - 




(Dl) 


The earth gravity field (per unit test mass) whose potential is Equation Dl is given as the gradient 
of Equation Dl, or 


r cos c p dk + r 


F 


(D2) 






or 



r 2 * 
+ 2 


+ 60(R 0 y/r) J 32 cos 2 sin 0 cos 2(\ “ \ 32 ) + 60 (R Q |r) J 33 cos 3 4 > cos 3(\ - \ 33 ) 

+ f (R 0 /r) 2 J 40 (35 sin 4 <f> - 30 sin 2 cf> + 3 ) 

25 

+ T ( R o/ r ) 2 J41 ( 7 s in 2 - 3 ) cos 4> sin <£ cos (X - \ 41 ) 

+ IT ( R o/ r ) 2 J42 ( 7 31112 ” 1 ) cos2 $ cos 2 (^ “ ^42) 

+ 525(R 0 /r) 2 J 43 cos 3 0 sin^>cos 3(^ - X 43 ) + 525(R 0 /r) 2 J 44 cos 4 <p cos 4(A. - X 44 ) 


(D 3 ) 


= 


Me 

r 2 


(R 0 /r) 2 | 6 J 22 cos<j!>sin 2 (^-^ 22 ) + 2 ( R o/ r ) J3 1 [ S sin2 ^ ~ *] sin ( X 'Si) 

3 °( R 0 / r ) J 32 cos<6sin^sin 2 (A.-\ 32 ) + 45(R 0 /r) J 33 cos 2 <$> sin 3(X - ^ 33 ) 

| (R 0 / r ) 2 J41 [7 sin 2 tf>- 3] sin 9b sin(\ - \ 41 ) + 15(R 0 /r) 2 J 42 (7 sin 2 4 > ~ l) cos 4 > sin 2 (\ - \ 42 ) 


+ 315(R 0 /r) 2 J 43 cos 2 4> sin sin 3(^ - A_ 43 ) 

+ 420(R 0/ /r) 2 J 44 cos 3 <t> sin 4(\ - ^ 44 )} 

~ (R 0 /r) 2 |-3J 20 sin<t> cos 4> + 6J 22 cos <p sin<£ cos 2^ - ^ 22 ) 

3 3 

- 2 ( R o/ r ) J 30 sln2 < t > ~ l) cos & + "2" ( R o/ r ) J 3 1 (15 sin 2 ~ ll) sin 4> cos (k - \ 31 ) 
+ 15(R 0 /r) J 32 (3 sin 2 l) cos <£cos 2(k~k 32 ) 

+ 45 (r o /r) J 33 cos 2 <£sin<£cos 3(\-X 33 ) - f (R 0 / r ) 2 J 4 o( 7 sin 2 </>-3) sin0cos<£ 

+ f- (R 0 / r ) 2 J 41 ( 28 sin4 0" 27 sin 2 <£ + 3 ) cos(\ - \ 41 ) 

+ 30(R 0 /r) 2 J 42 (7 sin 2 $ - 4) cos 0 sin 4> cos 2(k- \ 2 ) 

+ 105(R 0 /r) 2 J 43 (4 sin 2 l) cos 2 4>cos 3(k-k^ 3 ) 

+ 42O(R 0/ /r) 2 J 44 cos 3 <^sin^>cos 4(X.~ ^ 44 )} 


(D 4 ) 


(D 5 ) 
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II 


The actual sea-level surface of the earth is to be conceptualized through Equation D1 as a sphere 
of radius 6378 km, around which are superimposed the sum of the separate spherical harmonic 
deviations illustrated. To these static gravity deviations, of course, must be added a centrifugal 
earth- rotation potential deviation at the earth’s surface, to get the true sea level surface (see 
Reference D-l). 
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Appendix E 

Estimation of the Standard Error for a Function of Several Variables 
Whose Standard Errors, in an Experiment, are Known 


Let Z = z(x r X 2> • • • x n ) be the theoretical relationship between z and a number of variables 
X.. Let AX 1? ax 2 , • ■ *AX n be small errors of estimation in each of the X. variables. Then from 
the theorem for the total differential of a function of many variables: 


dZ 

az = 2^ ^r AXi • (E1) 

iTi 

Equation El represents the error made in Z for a single trial k at estimating the X i 's. Thus from 
one estimation k of the X i T s, Equation El gives 



i = 1 


Conceive now of a large number M of independent estimates k of the X i *s each giving a different 
error of estimation AX ik . One can ask from this extended experiment what will be the most likely 
value of the sums of the squares of the individual deviations AZ k . 

Thus, from Equation E2 



(E3) 
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II 


It is clear that if the AX ik T s are distributed with reasonable symmetry about zero, the value of 



AX ik AX jk 


will be close to zero for M sufficiently large. In fact, the expected or most likely value of this sum 
can probably be shown to be zero for any M when the AX ik f s are reasonably symmetrically distributed 
about the expected value of AX ik , which is zero. 

Using this result in Equation E3, it is found that 


E 


£ (AZk)2 



where E[ ] indicates the expected value of [ ], or 


E(AZ k 2 ) 


£$)'£>* 


) 


(E4) 


Equation E4 follows from the theorem that the expected value of a sum of random variables is the 
sum of the expected values of those random variables.* But the expected value of AZ k 2 is just the 
variance cr z 2 of z, and the E(AX i 2 k ) are just the variances of each variable X i? 

S («iJ = °X. ■ 


Thus Equation E4 becomes 


or 


Ma/ = 




(E5) 


* Bower, A. H., and Lieberman, G. J., Engineering Statistics , Englewood Cliffs, New Jersey: Prentice Hall, Inc., 1959. 
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The standard deviation of Z is thus 


a z = - 

In an experiment with a limited amount of data on the X. and <j x , if estimates of these, X. and S x , 
are obtained, then from Equation E6 an estimate of cr z is 


m 


(E6) 


S 


Z 



(E7) 


In the Syncom n data reduction for J 22 and k 22 by way of the ellipticity-drift parameters C,, C 2 , 
C 3 , the functional relationships between j 22 , k 22 and these parameters are (Equations 5, 6 and 2): 


k 


22 



J22 ^2 21 


A 22 = (c/ + c/)‘/> , 

- 72rr 2 (r q /a s ) 2 ' 2 


(cos 2 i s + l) 


(E8) 

(E9) 

(E9a) 


Applying Equation E7 to Equations E8 and E9 gives 


2(c/.C,>) [&>* K)’ * («.)“ (%)’]■ 

(c, 1 ‘C, 2 ) 122 [l C «)’ Kr * ( Sc ;)'j 


1/2 


(E10) 


(ElOa) 


Using the values of Sg^ , Sg^, c 2 and c 3 from the actual gravity experiment, Equations 40 through 
43, in Equations E10 and ElOa gives 


2 ( 


0.1 934 ‘,0.1 3^- ) [<°->«4 * 0.0481)> * (0. !355 * 0. 0145)’] 


1/2 


- 0-0852 radians - 4.87°, 


(Ell) 
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s* 


2 2 , actua 1 


10 4 (rad/day) 2 
(0. 19 34 2 + 0. 13 5 5 2 ) 1/2 


[(0.1934 xO-0145) 2 + (0.1355 x0.0481) 2 ] 


1/2 


- 3.00 x 10 6 (rad/solar day) 2 - 3.00 x 10 6 (rad/solar day) 2 x 0.9945 day 

= 2.98 x 10 -6 (rad/sid. day) 2 . 


Thus, 




/ [ . ( COs2 i s + ol 

S T 

J 22 , ac t u al 

( S A V 

\ A 22, actual// 

;727r 2 (R 0 /a s ) 2 1 2 J 


2. 98 x 10' 6 


13.859 [See (28)] 


= 0.215 x 10' 6 . 


(E12) 


The "best values of J 22 and \ 22 sensed from Syncom II drift over the Central Pacific with their 
standard errors given by Equations Ell and E12, are thus: 


J 22 = ( 1. 71 +0.22) x 10~ 6 

\ 22 = -(17.1+4.9)° . 


(E13) 
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Appendix F 

List of Symbols 


A 22 , j 22 , k 22 


a , 1 


c x . C 2 , c 3 


Second-order tesseral gravity constants of the earth 1 s field 

Semimajor axis and inclination of a 24-hour satellite orbit 

Constants in the drift equation for a 24-hour satellite, related to its mean 
motion and the two parameters of earth equatorial ellipticity 


R 0 Mean radius of the earth’s equator 

y Geographic longitude of the 24-hour satellite’s ascending equator crossing 
east of the minor equatorial axis of the earth 

cr2, s| f Z The variance and estimate of the variance of a random variable Z whose best 
estimate z is determined in an experiment 


NASA- Langley, 1966 
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